Extralimital terrestrials: A reassessment of range limits in Alaska’s land mammals

Understanding and mitigating the effects of anthropogenic climate change on species distributions requires the ability to track range shifts over time. This is particularly true for species occupying high-latitude regions, which are experiencing more extreme climate change than the rest of the world. In North America, the geographic ranges of many mammals reach their northernmost extent in Alaska, positioning this region at the leading edge of climate-induced distribution change. Over a decade has elapsed since the publication of the last spatial assessments of terrestrial mammals in the state. We compared public occurrence records against commonly referenced range maps to evaluate potential extralimital records and develop repeatable baseline range maps. We compared occurrence records from the Global Biodiversity Information Facility for 61 terrestrial mammal species native to mainland Alaska against a variety of range estimates (International Union for Conservation of Nature, Alaska Gap Analysis Project, and the published literature). We mapped extralimital records and calculated proportions of occurrences encompassed by range extents, measured mean direction and distance to prior range margins, evaluated predictive accuracy of published species models, and highlighted observations on federal lands in Alaska. Range comparisons identified 6,848 extralimital records for 39 of 61 (63.9%) terrestrial mainland Alaskan species. On average, 95.5% of Alaska Gap Analysis Project occurrence records and ranges were deemed accurate (i.e., > 90.0% correct) for 31 of 37 species, but overestimated extents for 13 species. The International Union for Conservation of Nature range maps encompassed 68.1% of occurrence records and were > 90% accurate for 17 of 39 species. Extralimital records represent either improved sampling and digitization or actual geographic range expansions. Here we provide new data-driven range maps, update standards for the archiving of museum-quality locational records and offer recommendations for mapping range changes for monitoring and conservation.


Introduction
Quantifying the degree to which species ranges have shifted over time is essential for understanding organismal responses to anthropogenic climate change.Although these efforts are key for biodiversity conservation [1], the process of documenting, quantifying, and predicting such shifts is challenging, particularly given the rapid pace of global change.For example, the effects of climate change at high latitudes are especially pronounced [2], with mean annual and seasonal temperatures increases of 2-4˚C over the past 20 years resulting in thawing permafrost, changing precipitation patterns, advancing shrub-and tree lines, and shorter fire-return intervals [3].The disproportionate effects of climate change at high latitudes are already causing a northward shift of the boreal forest at the expense of Arctic tundra, reallocating available habitat for Arctic-and boreal-adapted species [4][5][6].Although many species in these biomes are expected to maintain habitat associations and move vertically or laterally in accordance with shifting landcover [7][8][9][10], sufficient spatial data are lacking for many species [11], making range shifts challenging to document.Understanding the impacts of these habitat shifts on community assemblages requires accurate approaches for mapping and assessing species ranges.
Tracking species and community range shifts requires knowledge of prior species distributions and assessment of putative extralimital occurrences.Wildlife ranges are inherently dynamic, even during periods of relative stasis [12,13], so distinguishing true distributional shifts from newly collected information at poorly characterized range margins can be challenging.Thus, the identification of new extralimital occurrences is highly dependent on the range to which they are compared.One approach for developing a background range is to use species distribution models (SDMs) or ecological niche models (ENMs), which can incorporate environmental variables in often discontinuous predictions across landscapes [14][15][16].However, SDMs often lack repeatability, methods can vary considerably among studies and over time, and models may predict suitable but as-yet unoccupied marginal habitats [17,18].In contrast, the use of occurrences to delineate range extents (i.e., classic range maps) and identify extralimital records (i.e., occurrences outside identified range extents) may be the most robust and repeatable approach.Admittedly, range maps are coarse and generalized by nature and therefore often overestimate real occupancy [19], but they represent consistent baseline summaries of existing knowledge, often from peer-reviewed literature, against which new occurrences can be compared.
A further complication for identifying marginal shifts in species is that even verified occurrences indicate only that an individual of a species was found at a particular time and place and do not necessarily suggest regular habitation by that species [20].Furthermore, due to potential sampling biases or low marginal abundances, extralimital records could reflect improvements in sampling and knowledge of occurrences near range boundaries.Nevertheless, without developing accurate baseline range extents, monitoring species change across space will remain an unmet challenge.Ultimately, regular, geographically representative sampling that enhances comprehensive occurrence datasets is the best scientific means for monitoring changing range extents in the context of rapid climate change [1].
The ranges of numerous North American mammals reach their northern-and westernmost extents in Alaska, and so range limits in the state also take on continental and, for some, global significance [21].As we monitor changing species distributions and community assemblages, quantifying the dynamics of species near the leading edges of expanding ranges is paramount, as this space also represents the receding edge for several Arctic endemics whose ranges are likely contracting [4,[22][23][24].Climate change in Alaska is predicted to cause the geographic ranges of species to shift northward, upward in elevation, and towards the coast [4,23,25], yet the status of such predicted species-specific movements remains largely unknown.Limited research has examined the extent to which terrestrial mammalian range or distribution shifts have actually occurred for select species over historical timescales in Alaska, including the expansion of American Marten (Martes americana [Turton, 1806]) across the Kenai Peninsula [26], northward shifts in the Alaskan Hare's (Lepus othus Merriam, 1900) distribution in western Alaska [27], and the expansion of Snowshoe Hares (L.americanus Erxleben, 1777) [28], Moose (Alces alces [Linnaeus, 1758]) [29,30], and North American Beavers (Castor canadensis Kuhl, 1820) across the Alaskan Arctic [31].
Several research efforts conducted in the 2000s endeavored to collate and map occurrences of mammals in Alaska [21] and describe ranges of terrestrial vertebrates and other species (e.g., [32,33]).In their 2009 book, MacDonald and Cook [21] compiled information from georeferenced museum specimens for extant marine and terrestrial mammals in Alaska and mapped occurrences and marginal records across the state.This marked a notable milestone in Alaskan mammalogy, as no such effort had been attempted previously, and the resulting maps remain the standard against which new occurrences are interpreted.Between 2009 and 2014, the Alaska Gap Analysis Project (AKGAP) conducted a spatial review of distributions and range limits for all terrestrial vertebrates in Alaska using occurrence records, 4th level hydrologic units (HUs; represented by 8-digit hydrologic unit codes; HUC8), and three geospatial analyses [33].Independently, the International Union for the Conservation of Nature (IUCN) vastly expanded the scale and reach of its Red List accounts in the mid 2000's, digitizing and serving range maps at a global scale for nearly all species of mammals as well as more than 150,300 other species [32].
Since 2009, advances in digital data collection and the rise of online museum and citizenscience initiatives such as iDigBio and iNaturalist, respectively, have enhanced field detection of species and complemented the publicly available, georeferenced datasets from museums (e.g., Arctos, VertNet).Citizen-science applications such as iNaturalist allow observers to easily contribute georeferenced, photo-verified records of species from a variety of personal devices (complete with coordinates and metadata).If a plurality of users (regardless of their taxonomic expertise) agrees with the species identification, photo-observations may be deemed "research grade" and uploaded to the Global Biodiversity Information Facility (GBIF).GBIF is the world's largest and most comprehensive digital biodiversity repository (> 2.3 billion records) and it regularly harvests and collates records from iNaturalist and from dynamically curated, specimen-based databases like Arctos, VertNet, BISON (Biodiversity Information Serving Our Nation), and iDigBio, which have been steadily digitizing analog museum records.
The wealth of newly collected, newly georeferenced, and newly digitized specimens, as well as the growth in research-grade citizen science data, provides a timely opportunity to re-examine occurrence datasets of Alaskan mammals in relation to range estimates over the past decade.
Here we provide an update to the geographic status of northern terrestrial mammals with records beyond previously defined range extents.We review locational records [21,34] and range maps (IUCN and AKGAP) for 61 Alaskan mammals.We place new extralimital records in the context of federal lands; highlight species newly documented in Alaska; correct dubious records; evaluate the accuracy of commonly used range maps; contribute new, current, repeatable range maps; and review predictions of future range shifts made by previously developed distribution models.We also reassess the standards of identification, archiving, and public accessibility of georeferenced occurrence records.We undertook these analyses with the goal of improving the quality of source data underpinning range maps and predictive models so that occurrence record sets and current range maps can provide conservation planners and wildlife managers with improved means for monitoring changing terrestrial wildlife patterns.

Materials and methods
To obtain the most comprehensive set of publicly available mammal occurrence records, we searched GBIF on 25 February 2021 for all georeferenced (coordinate precision < 5 km) records (i.e., whole specimen, tissue sample, machine observation [remote detection], human observation) of terrestrial mammals known or suspected to occur in Alaska, USA [21,34].We quantified the total number of extralimital specimens for a species, as opposed to the number of extralimital localities, because the presence of multiple individuals at the same locality provides greater evidence of a robust extralimital population.We followed the taxonomy used by the American Society of Mammalogists Mammal Diversity Database (ASM; https://www.mammaldiversity.org;v1.11, released 15 April 2023 and accessed 23 February 2024) and evaluated 61 native, terrestrial, non-island-endemic mammal species (including bats) thought to occur naturally (i.e., not introduced) in Alaska (Tables 1 and S1; [21,33]).Baseline (2009) range polygons were constructed in ARCGIS 10.5 (ESRI, Inc., Redlands, CA) by digitizing marginal records for each species using coordinates provided by MacDonald and Cook [21] and then constructing minimum convex polygons (hereafter 2009 MCP) from the data.
In ARCGIS, we intersected species-specific GBIF record sets with their corresponding 2009 MCP to identify and tally extralimital records (those occurring at least 5 km outside baseline range margins; [21]) and used these species for further analysis (Table 1).We flagged records with suspected georeferencing or taxonomic errors as dubious and in need of further clarification, correction, or omission, and their institutions of origin were consulted accordingly.Records that were confirmed as accurate and those that were corrected by museum professionals met our criteria as "verified" and were included in final species occurrence datasets, whereas those remaining with implausible dates or locations were deemed unreliable and excluded from further analysis.
We focused subsequent analyses on those species with extralimital records (as defined in the previous paragraph; Table 1).For each species we intersected verified occurrences in ARCGIS with corresponding IUCN and AKGAP range maps [32,33] to determine the degree to which these maps corresponded to occurrence datasets.To calculate the sensitivity of IUCN and AKGAP ranges, we divided the number of occurrences within range extents by the total number of occurrences in Alaska (Table 1).For AKGAP ranges, we also calculated the percentage of HUs occupied within a range by dividing the number of HUs with at least one occurrence record by the total number of HUs in a range (Table 1).To remain consistent with prior research and to account for spatial errors related to coordinate precision, all intersections used a spatial tolerance of 5 km [21].We also identified species with new extralimital records that occurred in units managed by the U.We subsequently developed range maps for the 39 species with extralimital records, using our updated geodatabase of occurrence records.We used the Python tool, Concave Hull from Occurrence Points [35,36], to define concave polygons, representing current, inferred range extents for each species.Tool performance and range geometry is based on a minimum distance setting and a scale-dependent multiplicative factor that adjusts the aggregation distance to produce repeatable maps that do not over generalize [36,37].We used the default minimum search distance of 1,000 m, a buffer distance of 5 km, and a multiplicative factor of 12 for all species except Singing Voles (Microtus abbreviatus Miller, 1899) and Long-tailed Voles (M.longicaudus [Merriam, 1888]), for which we set the multiplicative factor to 8 and 15, respectively, to account for their disjunct distributions.

Range evaluations
From GBIF we harvested a total 124,425 records of 61 native, terrestrial, non-island-endemic mammals with occurrences in Alaska from 1818 to 2020 (Fig 1) and identified 39 species as having extralimital records beyond 2009 MCP range extents (Table 1) [21].Comparing occurrences for those species to IUCN and AKGAP range maps, we found that, on average, IUCN range maps encompassed 68.1% of GBIF records, whereas 95.5% of GBIF records were correctly included in AKGAP range extents (Table 1).AKGAP ranges exceeded 90.0% accuracy (correctly included presences / all occurrences) for 31 of 37 species (83.7%;AKGAP range maps were unavailable for 2 species of bats), whereas IUCN range maps reached 90% accuracy for just 17 of 39 species (43.6%) with available maps (Table 1).However, AKGAP ranges overestimated range extents for most species by including numerous HUs without occurrences.Among reviewed AKGAP range maps, an average of 57.3% of HUs contained GBIF records while HU occupancy exceeded 50.0% for 21 of 37 species (56.8%;Table 1).

Extralimital records
We identified 39 species as having �1 extralimital records since 2009 (Table 1 and Fig 2).From these we identified 6,853 extralimital records (including multiple specimens from the same locations) belonging to species from five orders (Eulipotyphla, Rodentia, Lagomorpha, Chiroptera, and Carnivora; Tables 1 and 2 and Fig 2 and S2    Recent extralimital specimens of Eurasian Least Shrews (Sorex minutissimus Zimmermann, 1780) and Eastern Pygmy Shrews (S. hoyi S.F.Baird, 1857) were collected in new areas across the Brooks Range (Noatak NP&P, Gates of the Arctic NP&P, and Arctic NWR [38]), exceeding their 2009 MCP ranges to the north by ~245 km and ~95 km, respectively.Barren Ground Shrews (S. ugyunak Anderson & Rand, 1945) were collected for the first time in Kobuk Valley NP and near Kanuti NWR over 130 km to the south of the 2009 MCP and were the only Arctic species in our analysis with recent occurrences along their southern range margin.
Recent extralimital specimens of other shrews were collected primarily in the vicinity of older records and so were not indicative of range changes.Extralimital specimens of Masked Shrews (S. cinereus Kerr, 1792) were collected in extreme Southeast Alaska, in Southwest Alaska as far as Unimak Island (Alaska Maritime NWR), on Nunivak Island (Yukon Delta NWR), and north of the Brooks Range (incl.Arctic NWR; Table 2 and Fig 2).Older extralimital specimens extended knowledge of their range to include the Beaufort Sea coast, as much as ~310 km to the north.Extralimital specimens of Tundra Shrews (S. tundrensis Merriam, 1900) were collected from the Seward Peninsula, in or near the Yukon Delta NWR, Alaska Peninsula NWR, Yukon-Charley Rivers National Preserve, and Arctic NWR, but also occurred near older specimens (Table 2 and Fig

Rodentia
We identified 5,446 extralimital records (including 1,312 after 2009) belonging to 17 rodent species that were distributed across the state (Tables 1 and 2     For some interior species with smaller ranges, recent extralimital records resulted in only small (< 10 km) extensions beyond the 2009 MCP range: Woodchucks (Marmota monax [Linnaeus, 1758]), Taiga Voles (Microtus xanthognathus [Leach, 1815]), and Northern Meadow Jumping Mice (Zapus hudsonius [Zimmermann, 1780]; Figs 2 and 4).Extralimital records in Southeast Alaska added new islands and areas of the mainland to the inferred ranges of rodent species whose extents terminate in this region, including Northwestern Jumping Mice (Zapus saltator Allen, 1899; Fig 2).While we identified 1,047 extralimital records of Northwestern Deermice, analyses for this species could not be separated from that of the closely related Eastern Deermouse given their fluid taxonomy.Genetic barcoding is necessary to accurately discern these species [39], and we had low confidence that GBIF records of P. maniculatus within the range of P. keeni were correctly identified.

Lagomorpha
We report 133 extralimital records for the three species of native lagomorph that occur in Alaska, including 55 since 2009 (Table 1).Extralimital records were located primarily around the periphery of the state ( Older extralimital specimens of Snowshoe Hares were collected across the North Slope (of the Brooks Range), extending their inferred range at least 75 km to the west.Several recent records were identified from the Seward Peninsula, Southwest Alaska, Kodiak NWR, and northern Southeast Alaska (Tables 1 and 2 and Fig 2).Numerous recent and older extralimital records of Alaskan Hares occurred in Northwest Alaska, including in Selawik NWR and Noatak National Preserve; [27]), resulting in a new range extent ~185 km to the northeast (Table 2 and Fig 2).Older extralimital records occurred at the base of the Alaska Peninsula (Katmai NP&P and Alaska Peninsula NWR), whereas newer records came from as far west as Unimak Island (Alaska Maritime NWR; Table 2).A single museum specimen (CMN 31930) believed to have been collected from near Utqiag ˙vik in 1898−well outside the Alaskan Hare's published or suspected range−was not included in our analyses for reasons discussed in Cason et al. [27].The Collared Pika (Ochotona collaris [Nelson, 1893]), whose range is limited to the southeastern quadrant of mainland Alaska, was represented by several recent extralimital specimens from the White Mountains, as much as ~45 km northwest of their 2009 MCP range, and one new observation (iNaturalist:7671018) from near Klondike Gold Rush NHP in 2017.Five specimens collected from Lake Clark NP&P in 2007 mark the westernmost known extent of this species in Alaska (Fig 2 ).

Chiroptera
Following the corrected identification of 15 Yuma Myotis (Myotis yumanensis [Allen, 1864]), specimens collected from southern Southeast Alaska [40] and the auditory detection of Northern Hoary Bats (Lasiurus cinereus [Palisot de Beauvois, 1796]) at four locations in northern Southeast Alaska in 2012-13, seven species of bat are now known to occur in Alaska (Table 1 and Fig 2).In total, we identified 159 extralimital records of bats (118 since 2009; Table 1).Except Little Brown Myotis (Myotis lucifugus [Le Conte, 1831]), all of these were from Southeast Alaska.On average, extralimital locations of bats were 95.  ).Little Brown Myotis is the only bat known to occur across Interior Alaska (between the Brooks and Alaska Ranges) and the Kodiak NWR [21,[41][42][43].The northernmost extralimital specimen (UAM: Mamm:127701) was collected in Wiseman, AK (adjacent to Gates of the Arctic NP&P) in 2013, pushing the extent of its known range ~170 km to the north compared to the 2009 MCP (Table 2 and Fig 2 and S2 Table).

Carnivora
We identified 436 records (74 since 2009) from six carnivore species as having extralimital records (Table 1 and Fig  American Black Bears (Ursus americanus Pallas, 1780), Red Foxes (Vulpes vulpes [Linnaeus, 1758]), and North American River Otters (Lontra canadensis [von Schreber, 1776]) had a mixture of recent and older extralimital records in both Southeast and northern Alaska (Table 1 and Fig 2 and S2 Table).Canadian Lynx (Lynx canadensis Kerr, 1792) and Red Fox were recently documented in the Alaska Peninsula NWR and Alaska Maritime NWR, respectively, demarking their southwestern range extents (Tables 1 and 2 and Fig 2).Three records of North American River Otters since 2014 in the Brooks Range and on the Beaufort Sea coast in the Arctic NWR are the northernmost records of this species in Alaska.Numerous extralimital records of Red Fox along the northern coast, especially 31 recent occurrences in northeastern Alaska near oil infrastructure and in the Arctic NWR, are representative of a range that extends 200 km to the northeast and now covers nearly all of Alaska.Numerous older extralimital records of Coyotes (Canis latrans Say, 1823) were identified from the Kenai NWR in the south but also recent records from Gates of the Arctic NP&P and Glacier Bay NP&P regions are indicative of (but likely still underrepresent) a Coyote range that has been expanding in Alaska for at least the past 70 years [21,44].Least Weasels (Mustela nivalis Linnaeus, 1766) occur across most of Western and northern Alaska, but a string of older specimens (including in Yukon Flats NWR) and new specimens in Eastern Alaska near Yukon-Charley Rivers National Preserve and Glacier Bay NP&P represent an eastward extension (~830 km) of the current inferred range of this species.

Range extensions
We identified extensions of previously inferred range limits in multiple directions for 39 terrestrial Alaskan mammals including two species (Northern Hoary Bats and Yuma Myotis) that were detected in Alaska for the first time during the past decade (Tables 1, 2 and S2).While new extralimital records may indicate the actual expansion of a range (e.g., Fig 3B), they could instead reflect previously unsampled localities where an understudied taxon may have long occurred.This can be a decidedly difficult distinction, especially given improved sampling in new and remote areas and the recent digitization of older records [45].This latter consideration best fits many of the historic specimen records not included in MacDonald and Cook's [21] review.During their assessment, many North American museums had either not yet digitized their collections or were not yet serving data to public online repositories.
After accounting for historically sampled yet previously unknown localities, anthropogenic climate change appears to be shifting the distribution of many species [8,46,47].Our primary aims here were to highlight extralimital records indicative of possible range extensions, to provide an update of the current geographic status of terrestrial mammals on NPS and FWS lands in Alaska, and to offer new, data-driven range maps that represent the current state of knowledge for terrestrial species at the forefront of climate change.While deciphering potential causes for any range shifts is complex and beyond our scope, we can compare observed extralimital occurrences against predictive geographic species distribution models (S1 Table) [22,25] and statistical habitat models [4] that predict where and how distributions of mammals are likely to change over the coming century.Below, we consider recent (i.e., after 2009) extralimital records in the context of available predictions of mainland terrestrial mammal species in Alaska for the 21 st century [4,22,25].

Eulipotyphla
Recent extralimital records of shrews were identified along Alaska's west coast and across the Brooks Range (Fig 2).Most extralimital specimens were collected in the vicinity of older records and likely represent more complete occurrence datasets rather than shifts in species ranges since 2009, especially because recent extralimital records indicate only modest differences in range extent versus the full extralimital dataset (Fig 3).However, specimens of Least, Pygmy, and Barren Ground Shrews were identified to the north and east of extents inferred from MacDonald and Cook's [21] maps.Hope et al.'s [25] models generally predicted nearstatewide (except the eastern Brooks Range and Northeast Alaska) distributions for Masked, Pygmy, Tundra, Least, and Western Water Shrews.Although models for these species were largely underpredictive for the current (2015) timeframe, model trends indicated expanding distributions for most shrews into at least the southern Brooks Range during the 2020s.Baltensperger and Huettmann's [22] models also predicted range expansions for these species onto the North Slope but contractions in Southwest Alaska by 2100.In contrast to Hope et al.'s [25] models, Baltensperger and Huettmann [22] indicate a small present-day (2015) range for Western Water Shrews that was predicted to expand northwestward across the state by 2100.Recent extralimital specimens for Western Water Shrews are consistent with this prediction.Both sets of models were in general agreement for most species, except Barren Ground Shrews.Baltensperger and Huettmann [22] predicted the range of this species to recede north of the Brooks Range by 2100, whereas Hope et al.'s [25] models predicted a gradual southward expansion of the range by the 2080s, a prediction consistent with the most recent extralimital record (MSB:Mamm:143102) for this species along the southern margins of the Brooks Range.

Rodentia
Observed extensions of previously inferred range margins (Figs 2 and 3) for the northwestern group of rodents are consistent with future distribution models that predict northwestwardly expanding distributions of Long-tailed Voles, Northern Bog Lemmings, Northern Red-backed Voles [22], and North American Red Squirrels [25].Statistical projections also predicted increased areas of available habitat for Northern Flying Squirrels, North American Red Squirrels, and Northern Red-backed Voles, whereas habitat areas for most other species were predicted to remain roughly the same [4].Hope et al.'s [25] models predicted North American Red Squirrels to occur on the Seward Peninsula by the 2050s, where an individual was recently observed near Nome, AK (although no other records have since been reported).Alaska Marmots were also recently detected in an area predicted by the 2080s [25] to be the last refugium for this species [25].Despite a predicted decline in habitat available for North American Beavers [4], several recent extralimital records and a comprehensive analysis of historical observations indicated that the species has, in fact, been expanding its distribution across northern and western Alaska [31].New extralimital specimens of Northern Bog Lemmings in Kobuk Valley NP were collected from areas where this species was predicted to occur over the coming decades [22,25].
Results for the second group of rodents (those with recent extralimital records primarily to the north and east; i.e., North American Porcupine, Meadow Vole, Singing Vole, Hoary Marmot), were mixed relative to prior model predictions.Newly documented occurrences of North American Porcupines 28−200 km north of their 2009 MCP range extent are consistent with Marcot et al.'s [4] projection that their habitat would increase over the coming century.Both Marcot et al. [4] and Baltensperger and Huettmann [22] predicted a decline in total available habitat and area of occurrence for Singing Voles and Meadow Voles.Despite overall predicted declines in Singing Vole distribution, Baltensperger and Huettmann's [22] models indicated a small northward extension of their 2100 distribution that includes the Beaufort Sea coastline, the same region from which several extralimital specimens were recently collected.Recent extralimital records of Meadow Voles ~25−112 km north of their 2009 MCP range lend support to Hope et al.'s [25] model that indicated an expanding distribution.Hope et al. [25] predicted a small northward expansion of Hoary Marmots by the 2080s, but models underrepresented occurrence records and did not include the area of eastern Interior Alaska where new extralimital records were recently documented ~35-135 km north of the 2009 MCP.

Lagomorpha
Extralimital records of Snowshoe Hares largely conformed to Hope et al.'s [25] distribution, though our current range map indicates broader occupation in Interior Alaska, a modest extension to the west and southeast (Fig 2).These results are consistent with trends of increasing shrub cover in tundra areas that should make these areas more hospitable for Snowshoe Hares [28], although this contradicts future habitat projections that indicated no significant change in available habitat [4].Recent extralimital records of Alaskan Hares indicated the occurrence of this species along the length of the Alaska Peninsula, and to north of the 2009 MCP [27].In contrast, Hope et al.'s [25] models did not predict the Alaska Peninsula to be part of the Alaskan Hare's range during the 21st century.New records from the Seward Peninsula are consistent with predictions of an expanding distribution across northern Alaska [27], but there were no recent extralimital records to indicate inland movement of their range as predicted by Hope et al. [25].In contrast, Marcot et al. [4] predicted a small and insignificant decline in available habitat for Alaskan Hares.Distribution models for Collared Pika predicted a constriction of their range over the coming decades [25], whereas recent extralimital records represent a small extension of their current observed range to the north and south, which may be a product of expanded sampling (Fig 2).

Chiroptera
Detection of bat species in Alaska has recently been augmented using acoustic recording to detect species in new locations, including the detection of Northern Hoary Bats, and the collection of Yuma Myotis in Alaska for the first time during the past decade.While new detections have expanded occurrence sets of bats in southeast Alaska [40,48], modeled trends in distribution, range, or available habitat at regional to statewide scales have not yet been conducted for mainland Alaska (but see [49,50]).This makes evaluating whether new extralimital records, especially to the north and west (Fig 3), represent actual range expansions or just improved sampling an even greater challenge.Given the small size of bat occurrence datasets, the difficulty in detecting bats in the field in many Alaskan habitats, and the technological developments that have made this work more feasible only recently, we suspect that new records of Chiroptera (more than other orders) represent a growth in the known occurrence set, rather than movements of species.More detection and modeling research is needed to improve our collective conception of the geographic status of bats in Alaska and we conjecture that further auditory detection efforts will continue to identify bat records in many new locations, especially on islands in Southeast Alaska.
Records harvested from GBIF generally corresponded well with MacDonald & Cook's [21] assessments, but new and newly digitized observations indicated possible range extensions > 5 km for 39 species.Over 1,700 recent extralimital records and their distance direction profiles (Fig 3) hint at the possibility that the range extensions summarized here may represent real movements and are not only artifacts of sampling effort or delayed digitization.Among initially identified extralimital records, just 4% contained errors and were flagged as dubious.However, we only checked records along range margins for accuracy, and numerous georeferencing and taxonomic identification errors likely remain in archives and online repositories.Additionally, given the lag between specimen collection, expert identification, and data digitization, new species occurrences or specimen-based range extensions resulting from recent field work and literature may not yet be reflected in recent GBIF data harvests.Finally, while museum and citizen science records provide a baseline for species presence, they are inapplicable to the documentation of species absences.Improving mechanisms for documenting and sharing absences of expected species are ongoing and would provide a better baseline for shifting distributional data.Only increased sampling and specimen archival from understudied regions of the state (e.g., Southwest and Western Alaska; Fig 1) will help to improve confidence that we are working with comprehensive and representative geospatial datasets, yet the difficult logistics and expense of conducting field work in remote Alaska have resulted in slow progress towards this goal.
To differentiate anecdotal observations from verifiable records, some type of physical record allowing for the positive identification of a species must be made available.We echo prior calls to explicitly recognize and rank evidentiary standards and advocate for the highest level of reliable evidence for the rarest species, while allowing other types of verifiable documentation for more common species [52].The most verifiable practice for identifying and archiving mammals is the collection and preservation of whole specimens, including skulls (dentition is necessary for differentiating many small mammals).This is particularly important for assessing species that have undergone taxonomic revisions or are difficult to identify, e.g., American Marten vs. Pacific Marten (M.caurina [Merriam, 1890]), Northwestern Deermouse vs. Eastern Deermouse.With the rise of digital citizen-science repositories, a sufficient level of verifiability may be met with a physical specimen, DNA, or eDNA sample, or quality photograph (although most North American small mammals cannot be reliably identified to species based on photos; [53]).For example, a photo that has been identified and confirmed by at least 2/3 of reviewers is necessary for iNaturalist records to become "research grade" and meet the standards for GBIF archival, although no formal training or expertise is required of iNaturalist reviewers.In some cases, diagnostic auditory evidence (e.g., for bats, canids, pikas, marmots, etc.) and photographs or impressions of tracks (e.g., carnivores, ungulates) are also sufficient and useful for the identification of some common or readily diagnosed species that may be difficult to observe otherwise.While some museums serve these types of observational data from their databases, community standards for distribution and access to disparate observational data types are still being developed.
We also advocate for adherence to the FAIR principles of scientific data management and stewardship: Findability, Accessibility, Interoperability, and Reuse of digital assets [54].To achieve a sustainable level of completeness for the global occurrence datasets of species, all natural history records should ideally be georeferenced, digitized, and served to a centralized public repository [55].Currently, over 1,600 institutions contribute records to one of these online repositories, but numerous federal and private institutions do not, at least on a consistent basis.To this end, iDigBio and BISON have been instrumental in digitizing older collections and datasets.Streamlining the dataflow from field observations to accurate, complete, and accessible digital collections is a goal that field biologists, management agencies, and natural history repositories alike can and should strive towards.Through this common workflow, we may improve confidence in the geospatial datasets that form the basis of range maps, distribution models, and quantitative analyses that inform management recommendations and conservation priorities.

Improving range map accuracy
Freely accessible range maps are a critical tool for assessing biodiversity patterns [56], but considerable gaps in coverage and questions about the accuracy of specific maps remain [37].For example, most range maps did not accurately reflect the extent of geographic records for 6 species (Keen's Myotis, Long-legged Myotis, Yuma Myotis, Long-tailed Vole, Northern Flying Squirrel, and Hoary Marmot; Table 1 and Fig 2).Further research into the recent biogeographic history of these species as well as new field-detection efforts are needed to further clarify and track changes.Ultimately, range maps both inform the public about species occurrence and are used extensively by scientists, who incorporate them into spatial analyses and use them to inform conservation assessments (e.g., [57]).As such, current and accurate data underpinning these maps is paramount.[21,[32][33][34].
Accurate range maps should attempt to balance errors of omission and commission.The AKGAP range maps [33] encompassed most species occurrences, but at the 4th level HU scale tended to overestimate occupancy.Regional GAP projects across the country have been valuable for identifying biodiversity hotspots and potential connectivity between disjunct distributions (e.g., [58]).While sensitivity of AKGAP range maps to occurrence data was generally quite good, occupancy % was poor (< 50.0%) for 13 species: Coyote, North American River Otter, Least Weasel, Little Brown Myotis, Collared Pika, North American Beaver, Singing Vole, Common Muskrat, North American Porcupine, Northern Flying Squirrel, Hoary Marmot, and Eurasian Least Shrew (Table 1).Such mismatches between occurrences and HU occupancy are the result of either a lack of records in under-sampled areas or range polygons that overestimated the geographic extent of these species [33].As such, we support the use of smaller hydrologic units (e.g., 6th level HUs; represented by 12-digit hydrologic unit codes; HUC12) to delineate range extents (e.g., [59]).In contrast, IUCN range maps encompassed fewer than 90.0% of occurrence records for a majority of analyzed species (Table 1).Methods guiding the development of IUCN range maps encourage species assessors to provide the best possible map based on available data [60], including occurrence point data, polygon data, and HUs [32].However, expert construction of range maps for many species did not appear to be well matched to the occurrence data in Alaska.Ensuring range maps used in spatial analyses, conservation assessments, and management decisions are accurate at statewide-to-regional scales has important relevance for the many downstream uses of these data.
Ideally, range maps will be updated in real or near-real time as new records are reported, so that archived maps may serve as a dynamic and reliable historic record of biogeographic change.This can be accomplished by improving species mapping and assessment processes so that species are regularly reviewed in a well-documented, referenced, and repeatable manner, which is key to providing better baseline biodiversity data.While having maps created or vetted by taxonomic experts can improve map quality [61], this is not always possible for all species.The option we implemented here was to analyze available occurrence records with a common range mapping tool with a limited number of settings [35].Consistent use of such a tool to map the concave range extents of multiple species while stipulating how parameters are set will help to remove subjectivity and overgeneralization from expert range maps and go a long way towards standardizing range maps across species and over time.

Conclusions
The actual and perceived ranges of species are continually changing, both as we assimilate occurrence records into increasingly networked and comprehensive data repositories and as species move in response to climate change and other anthropogenic and competitive stressors.Extralimital records and newly developed minimum concave range maps highlighted inferred range extensions as compared to 2009 MCP range polygons for 39 terrestrial mammal species in Alaska.As a disproportionate number of extralimital records were also identified from lands managed by the U.S. National Park Service and U.S. Fish & Wildlife Service in western and northern Alaska (especially Arctic NWR), these areas may serve as an early indication of changing mammal communities in Alaska as well as provide clear opportunities for monitoring species shifts.While it remains difficult to identify precise causes for changing species geographies, new extralimital observations were largely consistent with prior model predictions of northward and coastward movement of many boreal species over the coming century.To monitor changes in the geographic arrangement of species over time, it is imperative to continue field sampling efforts along range margins and to contribute voucher specimens and accurate occurrence data to appropriate public repositories so that we can use consistent range mapping methods (as proposed here) to regularly update species range maps and models over the coming century.
S. National Park Service (including National Parks [NP], National Preserves [which we do not abbreviate here for the purposes of disambiguation with National Parks], National Historic Parks (NHP), and National Parks and Preserves [NP&P]) and the U.S. Fish and Wildlife Service (including National Wildlife Refuges [NWRs]; Fig 1), two federal agencies whose holdings are uniquely situated to monitor range shifts on federal lands in Alaska.

Fig 1 .
Fig 1. Map of Alaska showing all records of native, terrestrial mammals downloaded from the Global Biodiversity Information Facility in relation to land units managed by the U.S. National Park Service and U.S. Fish and Wildlife Service (See Table2).

Fig 2 .
Fig 2. Maps of mammalian occurrence records.Maps for 39 species of terrestrial Alaskan species with extralimital records (those occurring > 5 km beyond marginal range polygons) harvested from the Global Biodiversity Information Facility (GBIF)[34] are included here (a-am) with range polygons from the International Union for the Conservation of Nature (IUCN)[32], the Alaska Gap Analysis Project (AKGAP)[33], and those based on the minimum bounding geometry of marginal records from MacDonald and Cook[21].Occurrence records are colored by year of collection with darker colors indicating more recent records and lighter colors indicating older records.Public lands managed by the U.S. National Park Service (green) and U.S. Fish and Wildlife Service (brown) are also shown.All maps are in the NAD83 Alaska Albers projection and have an extent of -173.304˚,51.139˚to -112.673˚,67.397˚, except maps of Southeast Alaska that are shown in reference to this extent (see inset).https://doi.org/10.1371/journal.pone.0294376.g002 2).Western Water Shrews (S. navigator [Baird, 1858]) collected in the White Mountains since 2006 added three extralimital specimens and increased their current range as much as 95 km to the northwest (Fig 2).
and Figs 2-4 and S2 Table).On average, all extralimital records of rodents occurred 47.9 km (SD = 38.2km) from 2009 MCP boundaries, at a mean bearing of 189.7˚(SD = 81.6˚;Figs 3 and 4), while extralimital records after 2009 occurred 50.9 km (SD = 31.0km) from 2009 MCP boundaries at a bearing of 210F ig 3. Wind rose of extralimital records.Wind rose graphs depicting the mean distance and direction of a) all extralimital records and b) recent (post-2009) extralimital records relative to each species' 2009 minimum convex polygon (MCP) range extent for the 5 taxonomic orders of terrestrial mammals considered (Carnivora = blue, Chiroptera = orange, Lagomorpha = purple, Rodentia = yellow, and Eulipotyphla = green).https://doi.org/10.1371/journal.pone.0294376.g003(SD = 62.3˚;Fig 3).While extralimital specimens indicated modest range extensions for rodents in several directions (Figs 2 and 3), different patterns emerged for two groups of species, one with range extensions to the west and northwest, and a second with extralimital records primarily to the north and northeast.

Fig 2 )
. Recent extralimital records of North American Porcupines extend the current known range ~200 km north of the 2009 MCP, across the Brooks Range, although still within the IUCN and AKGAP ranges (Fig 2) [32, 33].Thirty Singing Vole specimens collected from two locations near the Beaufort Sea coast since 2009 mark the northernmost records of this species, ~52 km north of their 2009 MCP range (Fig 2).Extensive sampling between 2010 and 2018 yielded 18 Hoary Marmot specimens from the Yukon-Tanana River Uplands and represent a ~100 km extension of their current range to the northeast (Fig 2).Similar collection efforts resulted in the detection of an Alaska Marmot (Marmota broweri Hall & Gilmore, 1934; UMObs:Mamm:146) in the Arctic NWR in 2011, extending the known range by ~30 km to the northeast (Fig 2 and S2 Table).NEON (National Ecological Observation Network) surveys also detected the northernmost occurrences of Meadow Voles at the University of Alaska's Toolik Field Station and in Arctic NWR in 2019, a northerly extension of their current range by ~105 km (Fig 2).

Table 1 .
(Continued) [35]ies attributes include the American Society of Mammalogists (ASM) taxonomy, the number of occurrence records, extralimital records (with those since 2009 in parentheses), range map sensitivities, specificities, and current (2020) range areas.Extralimital records are defined as those > 5 km outside of minimum bounding geometry of marginal records from MacDonald and Cook[21].We report calculations of the percentages of GBIF points located within IUCN and AKGAP range polygons as well as the percentage of 4th level hydrologic units (HUs) in AKGAP ranges containing at least one GBIF record.2020 Concave Range Areas (km 2 ) were calculated using a minimum concave range mapping tool[35]; * indicates ranges with different multiplicative factors.† indicates a single extralimital record of Odocoileus hemionus hemionus, a subspecies of the Mule Deer (O.hemionus [Rafinesque 1817]), in Interior Alaska, independent of O. h.sitkensis (Sitka Black-tailed Deer) populations elsewhere in Alaska.https://doi.org/10.1371/journal.pone.0294376.t001

Table ) .
Among these, 1,702 extralimital detections occurred in 2009 or later.Of the extralimital records, 3,932 belonged to a single species, the Northwestern Deermouse (Peromyscus keeni [Rhoades, 1894]), whose taxonomy is entangled with that of the Eastern Deermouse (P.maniculatus[Wagner, 1845]).We flagged an additional 310 records (4.5%) that appeared to be extralimital, but which were incorrectly georeferenced or misidentified taxonomically.Through communication with museums, we resolved georeferencing errors for these specimens.Below, we describe the extralimital dataset by taxonomic order with a focus on how recent (� 2009) extralimital records compare to prior range limits (Table1 and Fig 2) and how they are distributed across lands managed by the U.S. National Park Service and U.S. National Wildlife Service (Table2and Fig1).Additional information and datasets are available in the Supplementary Files (S1 and S2 Tables and at Open Science Framework (https://osf.io/h6u5n).
EulipotyphlaA total of 677 shrew records, belonging to six species, were identified as extralimital (occurring outside 2009 MCP species ranges), including 141 that were collected in 2009 or later.On average, extralimital specimens of eulipotyphlan species occurred 48.0 km (SD = 71.9km)from their nearest 2009 MCP boundary at a bearing of 236.5˚(SD = 100.1˚;Fig3A).Recent extralimital shrews were collected an average of 36.2 km (SD = 33.1)from the nearest 2009 MCP boundary of each species at a bearing of 185.0˚(SD = 120˚; Fig 3B).Recent extralimital records primarily occurred across the northern half of the state, but the most distal extralimital specimens were identified in Southeast and Southwest Alaska before 2009 (Table 1 and Figs 2-4 and S2 Table).

Table 2 )
. Northern Flying Squirrels, commonly associated with eastern Alaska, were collected as far west as McGrath, AK in 2013 (UAM: Mamm:117637; see [21] for abbreviations used for museum collections) and west of Denali NP&P between 2014 and 2016, extending the known range of this boreal species ~120 km westward (Fig 2 and S2 Table).A North American Red Squirrel (UAMObs:Mamm:254) was photographed on the western Seward Peninsula in 2019, marking the westernmost extent of this species, normally associated with the boreal forest (S2 Table).Two Northern Bog Lemming specimens (MSB:Mamm:291027, 291034) collected from Kobuk Valley NP in 2014 represent the northwesternmost detections of this species and extend the current known range extent by ~160 km (Fig 2 and S2 Table).Extralimital records for a second group of rodents extended inferred ranges primarily to the north: North American Porcupine (Erethizon dorsatum [Linnaeus, 1758]), Singing Voles, Hoary Marmots (Marmota caligata [Eschscholtz, 1829]), and Eastern Meadow Voles (Microtus pennsylvanicus [Ord, 1815];